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Abstract 

We consider the question of linear instability of an equilibrium of the 
Relativistic Vlasov-Maxwell (RVM) System that has a strong magnetic 
field. Standard instability results deal with systems where there are fewer 
particles with higher energies. In this paper we extend those results to 
the class of equilibria for which the number of particles does not depend 
monotonically on the energy. Without the standard sign assumptions, the 
analysis becomes significantly more involved. 



1 Introduction 

1.1 The Relativistic Vlasov-Maxwell System 

In this paper we consider the linear stability of a super-heated plasma. The 
plasma is assumed to have low density, and thus collisions between particles 
may be ignored. The behavior of such a system is governed by the Relativistic 
Vlasov-Maxwell (RVM) System of Equations 

9 t / ± + v • V x /± ± (E + W xt + v x (B + B ext )) ■ V v f ± = 0, (1.1a) 
d t E = VxB-j, V-E = p, p = J(f + -f-)dv, (1.1b) 

9tB = -VxE, V-B = 0, j = j «(/+-/-) dv, (1.1c) 

where v is the momentum, v = vj yl-j- |t>| 2 is the velocity and all physical 
constants (like the speed of light c and the mass of the particles m) arc set to 
be 1. The superscripts ± indicate the two species - electrons and ions. The 
transport equation (|l.la|) is called the Vlasov Equation, and it is coupled with 
Maxwell's Equations. E(t, x) and B(t, x) are the electric and magnetic fields, 
/ ± (f, x, v) > are the electron and ion distribution functions, and E ext and 
TS ext are external electric and magnetic fields. In addition, p(x) is the net 
charge density at the point x, and j(x) is the current density at the point x. 
Particle charges are always +1 or — 1. The assumption that all masses are equal 



is physically relevant: It often occurs that interstellar plasma is made up of dust 
particles of uniform mass charged positively and negatively. 

The subject of this paper is the "purely magnetic" case, where the electric 
field is assumed to be negligible compared to the magnetic field, at least on the 
typical scale of the system. Together with a certain symmetry assumption on 
the two species, this assumption makes certain aspects of the analysis simpler, 
while keeping most of the interesting physical aspects. This result is outlined 
in Section 9 of [33], though in [33] the authors fail to consider all possible cases, 
and, therefore, their proof outline is incomplete. 

Our equilibrium is represented in terms of two invariants of the particles - 
the energy e and the momentum p. We explore the question of linear instability, 
and give concrete criteria for instability in terms of the spectral properties of 
certain Schrodinger operators acting on periodic functions of x alone, not v. 

1.2 The l| Dimensional Case 

For simplicity, we take the lowest dimensional system which has a nontrivial 
magnetic field and which has many physical applications: The so-called \\ 
dimensional case. In this setting, we have one spatial dimension and a two- 
dimensional momentum space. The single spatial variable x corresponds to v\, 
and the additional velocity dimension is denoted by V2 ■ We write v = (v\ , vi ) . 
It turns out that this dimensional setting is actually quite significant even for 
physicists: When one considers a tokamak (which typically has cylindrical sym- 
metry), then locally, in a small region, the symmetries that one observes are 
precisely those described here. 

The notation in the \\ dimensional case is as follows: We let f ± (t,x,v) 
be the ion and electron distribution functions, E(t,x) = (Ei(t, x), E-zit, x), 0) 
and ~B(t,x) = (0, 0, B(t, x)) be the electric and magnetic fields. In addition, we 
define the electric and magnetic potentials <f> and ip, which satisfy 



We assume the existence of an equilibrium f 0,± (x,v) which is a solution of 
RVM. By Jeans' Theorem (cf [6]) this equilibrium can be represented in the 
coordinates 



as / 0,± (a;, v) = /i ± (e ± ,p ± ), where (fP and ip° are the electric and magnetic 
potentials, satisfying E® = —d x (f> and B° = d x ip°, with E® and B° being the 
equilibrium electric and magnetic fields. In addition E% = 0. We note that 
the purely magnetic assumption implies that tfi° = 0, and, hence, e + = e~ 
can simply be denoted by e. Henceforth it will be understood that ^ are 
evaluated at e 1 p ± , so we will simply write ^ ± (e,p). In this paper, we consider 
the "nonmonotone" case; by that we simply mean that 



d x <t> 
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on some subset of the set {^{e.p) > 0}. Here, "fif" means "the derivative of 
/i* with respect to the first component evaluated at (e 7 p ± ). v Similarly, "/Zp " is 
the derivative with respect to the second component. By "monotone" , we mean 
that fif < on the set {p^ > 0}. Roughly speaking, the coordinates e and p 
should be understood to be energy and momentum respectively. 



1.3 Main Results 

We assume for simplicity that the equilibrium has some given period P. We 
define two operators acting on functions of x, whose properties will be rigorously 
treated later, as 

A\h = -d 2 x h - J W dvj h + J2 J vt^h dv 

A° 2 h = -d 2 x h dv ^j h -J2j tfW^h) dv, 

where V^- are projection operators onto some subspaces of a certain Hilbert 
space, also rigorously defined later. The operators A\ and A% have as their 
domain the space 

H 2 P = {his P periodic on R and h e H 2 (0, P)} 

where H 2 (0, P) is the usual Sobolev space of functions on (0, P) whose first two 
derivatives are square integrable. We also define the number 1° as: 



(vi) dv dx. 



In order to properly define function spaces that include functions that do not 
necessarily decay at infinity, we define a weight w that has the form 

w = c(l + \e\y a (1.4) 

for some a > 2 and c > 0, where e is the particle energy as defined above. We 
require /i ± <G C 1 to satisfy 

(\l4\ + \l4\)(e t p)<w(e) (1.5) 

so that f(\nf\ + \fif |) dv is finite. Obviously, we require // ± (e,p) > 0. 

As previously mentioned, we assume that the equilibrium is purely magnetic 
- that is, we assume that it is arranged in such a way so that there is no electric 
field, only a magnetic field - and satisfies a symmetry property that ensures 
that p° = 0: 

0° = O and fj, + (e,p) = fj,~(e,-p) (1.6) 

Lin and Strauss prove the existence of such equilibria in the appendix of [34] . 
The assumption <fi° = is physical: In many physical systems the equilibrium 
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state has negligible electric fields which are only significant on small scales, while 
the magnetic field is significant on the typical scale of the system. 
In this paper, we denote 

neg(J-) — {the number of negative eigenvalues (counting multiplicity) of the operator J 7 } . 

Similarly, pos^J 7 ) denotes the number of positive eigenvalues. Later we will 
show that A\ and A® have discrete spectra. For a real number a, neg(a) is 1 if 
a < and otherwise. 

Theorem 1. Let f°'^(x, v) = /i ± (e,p) be a periodic equilibrium satisfying (|1.5[) 
and (|1.6[) . Assume that the null space of A® consists of the constant functions 
and that I ^ 0. Then the equilibrium is spectrally unstable if 



Theorem 2. Under the additional assumption that the null space of A\ is 
trivial, the equilibrium is spectrally unstable if 



By "spectrally unstable" we mean that the system linearized around the 
equilibrium solution /i ± (e,p) has a purely growing mode solution of the form 



The first step in proving these results is linearizing the problem (see §2.2p . 
We then proceed by integrating Maxwell's Equations along the particle trajec- 
tories. The immediate benefit of this, is that we lose all dependence on the 
variable u, and can work, for the remainder of the paper, with functions of 
x alone. These calculations appear in £12.41 where operators that are closely 
related to A\ and A\ naturally arise. 

The integration along the particle paths leads us to three equations in three 
unknowns. Moreover, these equations depend upon the exponential growth pa- 
rameter A: finding a growing mode is equivalent to finding a nontrivial solution 
to this system for some < A < oo. For convenience, we write the system as a 
matrix of operators, denoted by M. x . We show that this matrix of operators is 
selfadjoint. Our general method for showing that the system has a solution for 
some < A < oo is to keep track of the spectrum of M. x as A varies, and show 
that it must cross 0. 

However, the matrix operator Ai* has point spectrum extending both to — oo 
and to +oo, due to the appearance of both d% and —d%- Therefore, to keep track 
of the spectrum, in Sj3]wc truncate the matrix operator into a finite-dimensional 
mapping that we denote by where n is the truncation parameter. A first 
obstacle is to show that for fixed n the spectrum of M.^ varies in such a way 
that a nontrivial kernel is guaranteed for some < A„ < oo. We do this by 
showing that has a different number of negative eigenvalues for small values 
of A and for large values of A and then applying the intermediate value theorem. 



neg (A? 2 ) > neg (A°) + neg(-l°). 



neg (A^)^ neg (A$) + neg(-l°). 




A > 0. 



(1.7) 
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The next difficulty lies in showing that the A n 's are uniformly bounded away 
from and from +00 for all n. We show this in Proposition [376] and in Lemma 
1531 

In ^~4]we recover our original problem, by "taking the limit n — ¥ 00." This 
process is tricky in itself, since we must define a proper sense in which finite- 
dimensional mappings "tend" to an infinite-dimensional mapping on Hilbert 
spaces. This is discussed, in part, in Lemma 13.41 appearing in the preceding 
section. 

Once we have shown that A4 X has a nontrivial kernel for some < A < 
00, in SjE] we demonstrate that the resulting solution is indeed a solution to 
the linearized RVM system. SJ7] contains many of the technical details used 
throughout the paper, whereas in fj6]we demonstrate our results in two examples. 

The main difficulty of this paper is in the truncation, and in recovering the 
original problem from the truncated one. Moreover, the truncation is performed 
by means of a spectral projection onto finitely many eigenvalues of How- 
ever, the spectrum of Ai x may change quite dramatically at A = 0, since all 
operators only converge strongly as A — > 0, but not in operator norm. There- 
fore, an additional difficulty is the need to deal with two varying parameters 
simultaneously - the growth parameter A and the dimension of the truncation 
n. 

Our result extends the results of [33], and, more specifically, [53], where 
Lin and Strauss established a linear stability criterion for equilibrium solutions 
that arc strictly monotone, that is, the number of particles at a given energy 
strictly decreases with the energy. They define a selfadjoint operator £°, acting 
on functions of x, with a rather simple spectral structure, and derive stability 
conditions based on those spectral properties. Namely, they prove that C° > 
implies spectral stability (in [53]), while C° ^ implies that there exists a grow- 
ing mode (in [33]). The monotonicity assumption, which we drop in this paper, 
proves to be important in [34] [33], as it determines the signs of some of the 
operators. This implies that £° is not relevant to our case. The root cause of 
most of the analytical difficulties of our paper is the inability to determine the 
signs of any of the operators. 

This paper is organized as follows: 

1. Introduction. 

2. Setup. We present our function spaces and operators. 

3. Behavior For Small and Large A. We "count" eigenvalues for small 
and large values of A. 

4. Limit as n — > 00. We retrieve our original system from the finite- 
dimensional approximations. 

5. Construction of a Growing Mode. We verify that indeed the system 
retrieved is the right one. 

6. Examples. We give two examples that illustrate our results. 

7. The Operators. We discuss the main properties of the operators. 
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2 Setup 



2.1 The Function Spaces 

The function spaces we use throughout this paper are as follows: We denote 
by I/p and Hp the usual Sobolev spaces of P-pcriodic functions on R, that are 
square integrablc on the interval [0,P], as well as their first k derivatives. In 
addition, wc denote: 



T 2 



H 



PA) 



T 2 




m{x, v) m is P-periodic in a;, ||m| 



|m| 2 w dv dx < oo 



where w is the weight defined in (|1.5[) . We note that L 2 W is independent of the 
choice of sign ±, since w only depends upon e which is independent of ±. In 
addition, we use the following notation: 



1^2 and (■, ■) L 2 denote the norm and inner product in L 2 P respectively 
■\\ w and (•, •) denote the norm and inner product in L 2 ^ respectively 



2.2 The Basic Equations 

In the 1 h dimensional case the RVM system becomes a system of scalar equa- 
tions: 

dtf± + vAf^ ± {E 1 + t) 2 S)^ 1 / ± ± (P 2 - t' 1 B)9„ 2 / ± = (2.1a) 

dtEi = -j x (2.1b) 

d t E 2 + d x B = -j 2 (2.1c) 

d t B = -d x E 2 (2.1d) 

d a .E 1 =n + p (2.1c) 



where v = vj (v) and (v) = yl + \v\ 2 , and 

P = j (/+ - r ) dv (2.2a) 
ji= f Hi (/+-/") dv, i = l,2 (2.2b) 
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and the external fields are replaced by the constant external radiation hq . From 
our definitions (|1.2[) of the electric and magnetic potentials, we get the second 
order ODEs for <f) and tp'- 

%</> = -d x Bk = -p=- j " (/+ - /") dv (2.3) 

d 2 x i> = d x B = -j 2 - d t E 2 = -[v 2 (f+- /") dv - d t E 2 . (2.4) 
The linearized Vlasov equation is 

(d t + D±) /± = ^nt<kEi ± pj^ifl T (^«a + ) #2, (2-5) 

where 

D ± = vA ± (E° + v 2 B°) d Vl T ViB°d V2 . 

It is important to note that the purely magnetic assumption implies that 
E® = 0, but the perturbation E\ need not vanish. Since we are looking for a 
(purely) growing mode with exponent A (see l|1.7|) ). we replace everywhere dt 
by A. Thus, our equations for the P-periodic electric and magnetic (perturbed) 
potentials (f> and ip become 

B = d x i\> 

along with 

E 2 = -\ij>, 

which is a result of the integration of (|2.1d[) and setting the constant of inte- 
gration to be 0, and, finally, 

Ei = -d x <f> - Xb, 

where b € R is meant to allow E\ to have a nonzero mean. Replacing d t by A, 
the linearized Vlasov equation becomes 

(A + D±) /± = ±fxfv 1 (d x <l> + Xb) ± (i^vAi) ± A (fifv 2 + 4) ^ ( 2 - 6 ) 
along with Maxwell's Equations 

~ Xd x ct) - X 2 b = XE X =-ji (2.7a) 

- X 2 ip + d%<ip = XE 2 + d x B = -j 2 (2.7b) 

- 8l<j> = d x E\ = p. (2.7c) 

We see that in these equations there is only dependence upon derivatives 
of the electric potential 0, and never dependence upon <j) itself. Therefore, 
throughout this paper we let (j) G L 2 P0 . Now we introduce the particle paths 
(X ± (s- 1 x, v), V ± (s; x, v)) of the equilibrium state, governed by the transport 
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operators D ± , where — oo < s < oo. They satisfy the system of ordinary 
differential equations 

X ± = V^ (2.8a) 
V± = ±V 2 ± B°(X ± ) (2.8b) 
Vf = tV^B ^), (2.8c) 

where □ = d/ ds is the derivative along the characteristic curves, and the initial 
conditions are 

(X ± (0,x,v),V ± (0,x,v)) = (x,v). (2.9) 

When there is no risk of confusion, we abbreviate X ± = X ± (s) = X ± (s;x,v) 
and = V ± (s) = (Vf 1 (s; x, v), V 2 (s; x, v)). Now we rewrite the Vlasov Equa- 
tion integrated along the particle paths. Here it is crucial that e^ 1 and p^, and 
any function of these variables, are constant along the trajectories. This implies 
that, /Zg and /i^ are constants under s-diffcrcntiation. We multiply (|2.6p by 
e As and notice that the left hand side becomes the perfect derivative (e As /). 
Integrating the right hand side along the particle paths, one has 

J — OO 

l4 (vtd^X^) + A0(X ± )) ds t f Xe^nfftX*) ds 

^ ^ J — OO 

± / e Xs ^ (yfd^ix*) + Mix ± ) \ ds 

J — oo 

± f Xe Xs nf (vfb + V 2 ± ip(X ± )') ds = I + II + III + IV. 

J — oo 

Recalling that D + and D~ both reduce to v\d x when applied to functions 
of x alone (and not v), we see that the integrands in terms / and III are 
e Xs /j,f (D(f) + X(j>) and e Xs /j,f (Dip + Xip), respectively, evaluated along the ap- 
propriate particle paths. Therefore, both integrands become no more than 



= ± / e**-* 



[e" s fj,f(j)) and ^ (e Xs fif r ip). We conclude that the terms I and III become 



±e //g 4>(x) and ±e [Apip(x), with no boundary terms due to our decay as- 
sumptions. The other terms arc kept in integral form as above. Since /«;f are 
constant along the trajectories, we may evaluate them at s = 0, and they have 
no role in the integration. After dividing both sides by the exponent, we finally 
have 

f ± (x,v)=±£<Kx)± f j±1>(x) (2.10) 



ds. 



We simplify this expression by introducing the operators Q± : L 2 W — > L^, 



defined by: 

,o 

(Q±k)(x,v)= / Xe Xs k(X ± (s;x,v),V ± (s;x,v)) ds 
where k = k(x, v) : [0, P) x R 2 ->• R. 

Remark 2.1. We note that if h{x) £ Lp, then h(x,v) := h(x), (x,v) S 
[0, P] x R 2 , is clearly in L^. Therefore, Q± act on functions in L 2 P as well. 

With our definition of Q± , (|2.10j) becomes 

f + (x,v) = +v+(/>(x)+^(x)-v+Q*(p + ii+Q^(hi>) + b(itQ+vi 

f~(x,v) = - / j l -(l)(x)- f i-i;(x)+ f i-Q x <p-fx-Q x (v 2 '>P)-b f j,-QH 1 . 

(2.11) 

An important quantity which we use frequently, is / + — / . We write it 
explicitly for future reference: 

/+ - r = £ 0** + - ^Q±^ + £ol(w) + b^Qih) 

(2.12) 

We make the following two useful observations: 

1. It holds that 

(j£ + v 2 vt) dv = J^- dv = 0. (2.13) 

2. The particle paths preserve volume. That is, for fixed s, 

(X (s; x,v),V (s; x,v)) — > (x, v) both have Jacobian = 1. (2-14) 

Lemma 2.2 (Properties of D ± ). D are skew-adjoint operators on L^. Their 
null spaces ker D ± consist of all functions g = g(x,v) in that are constant 
on each connected component in R x R 2 of {e = const and p = const}. In 
particular, ker D ± contain all functions of e and of p . 

Proof. We show for the '+' case, and drop the + superscripts. It is straightfor- 
ward to verify that De = Dp = 0. Therefore kerD contains all functions of e 
and of p. Skew-adjointness is easily seen due to integration by parts, as D is a 
first-order differential operator. Derivatives that "hit" w vanish, since w = w(e) 
is a function of e. □ 

Definition 2.3. We define V ± to be the orthogonal projection operators of I? w 
onto kerD*. 

Lemma 2.4. The projection operators V preserve parity with respect to the 
variable v\. 
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Proof. Let us demonstrate for V + and drop the + superscript to simplify nota- 
tion. The demonstration for V~ is identical. Recall that 

D = D+ = vid x + v 2 B°d Vl - viB°d V2 . 

Now, let / = f(x, Vi, V2) and let R be the operator that reverses V\ : Rh(x, Vi, V2) = 
h(x, — V\, ua). Then 

D(Rf) = -R{vid x f) - R(v 2 B°d v J) + Riv^dvJ) = -R(Df) 

Therefore / £ ker D if and only if Rf £ kerD. This implies that one can find 
a basis of even and odd functions (in the variable t>i) to the space kerD. To 
show that if / is even (odd) in V\ then Vf is also even (odd) in vi, we let 
g £ ker D be, without loss of generality, even or odd. Then it must hold that 
//(/ — V f)g w dx dv = 0. In the case that / is even, we change variables 
vi -> -Vi to get ±JJ(f - R(Vf))g w dx dv = and, therefore, R(Vf) = Vf. 
Here the ± depends on the parity of g. The odd case is treated in the same 
way. □ 

Lemma 2.5 (Properties of Q±). Let < A < 00. 

1. Q± map L 2 W — > L 2 W with operator norm = 1. Moreover, recalling Remark 
12.11 Q± are also bounded as operators L? p — > L 2 W . 

2. For all m(x,v) £ L 2 W , \\ Q^m - V ± m\\ w -> as A -> 0. 

3. If a > 0, then \\Q\ - Q±\\ = 0(|A - cr\) as A -> a, where || • || is the 
operator norm from to L^. 

4. For all m(x, v) £ L"^, \\Q±m — m|| — > as A — > 00. 

5. Let v = (—vi, V2) and let fi{x, v) = n(x, v). Then (Q^.m, n) = (m, Q±/i) . 

Proof. Let us show for the '+' case, and drop the + indices to simplify notation. 

1. We demonstrate showing the two bounds on Q x using a dual method and 
using a direct method; this shows that either method works, and each has 
its benefits. Let m,n £ L 2 W . 

(Q X m,n) w = [ \e Xs [ f ' (my/w) (X(s),V(s))-(jiVw)(x,v) dv dx ds < \\m\\ w \\n\ 



—00 



where in the equality we used the fact that w — w(e) is constant along the 
trajectories (this is true since e is constant along the trajectories), and the 
inequality is simply the Cauchy-Schwarz Inequality. The assertion that 
the operator norm is 1 is verified since Q x l = 1 (note that 1 £ L"^). 
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Now, for the L 2 P bound we let h(x) £ L 2 P and use (|2.14p to get 



\Q h\ 



Xe As h(X(s;x,v)) ds 



< 



— oo 




if dv dx 



Ae s h(x)\ ds I w dv dx 



(2.15) 



w dv \h(x)\ 2 dx 



< ( sup I w dv) \\h\\ 2 L 2 = C \\h\\ 2 L 2 



2. We let M denote the spectral measure of the selfadjoint operator T = 
—iD in the space L 2 W . We can therefore write, for any m = m(x,v), 
m(X(s), V(s)) = e sD m = e tsT m, so that 



Q m 



\e* s m(X(s),V(s)) ds = / Ae / e tas dM{a)mds 



X 



A + ia 



dM(a)m. 



The projection V can be written as V = M({0}) = J R X dM where x(0) 
1 and x(a) = for a £ 0. Thus 



\Q x m-Vm\\ = 



Xe Xs m(X(s),V(s)) ds~Vm 



X 



X + ia 



d\\M{a)m\\ 



due to the orthogonality of spectral projections. We use the dominated 
convergence theorem to finish the proof that, indeed, this expression tends 
to as A -> 0. 

3. To show that \\Q± — Q±\\ = 0(\X — a\) as A — > a we use the fact that 
(x,v) — > (X, V) has Jacobian= 1: 



\Q m — Q a Tn\\ w < 







<- f 


Ae As - ae as 


J — oo 








- /" 


Xe Xs -ae as \ 


J — oo 





< C|lnA-lno-| \\m\\ w . 
4. With M the same spectral measure as above, we have: 



Therefore 



Q m — m 



\Q x m-m\\l< 



X 



X + ia J 



1 ) dM(a)m. 



X 



1 



A + ia 

as A — > oo, by the dominated convergence theorem. 



d\\M(a)m\\l 
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5. 

(Q x m,n) w = Q x (m(x,v)) n(x,v) w dv dx 



o 

J (^J \e Xs m(X(s), V(s)) ds^j n(x, v) w dv dx 
p , ,o 

/ / Xe Xs m(x, v) n(X(—s),V(—s)) w ds dv dx 

J J-oc 

P r / r0 \ 

\e Xs h(X(s), V(s)) ds I m(x, v) w dv dx 

— oo / 



where for the third equality we used (|2.14[) and the fact that w is invariant 
under D, and for the fourth equality we used the fact that 

X(-s;x, -v 1} v 2 ) = X(s;x,v 1 ,v 2 ) 
-Vi(-s]x,-V!,v 2 ) = Vi(s;x,v 1 ,v 2 ) (2-16) 
V 2 (-s;x, -Vi,v 2 ) = V 2 (s;x,vi,v 2 ). 

This type of calculation appears often, so it is worthwhile writing in detail. 

□ 

2.3 The Operators 

In addition to the definitions of A® and A^ in ^1 .31 we define the following 
operators depending on a real parameter < A < oo, acting on L 2 p^,L 2 P ,L 2 P , 
with domains Hp ,Hp,L P , respectively: 

A x h = -%h - (E / £ du ) h + E / £ Q ± h dv > 

A x h = -d 2 x h + X 2 h -{jlj dv j h - E / ^ v 2 Q 

B X h = (E / ^ d ^j k + T.j VeQii^h) dv. 

As we have seen in Lemma |2~51 — > V ± strongly as A — > 0. We also define 
the following multiplication operators with domain K and range L 2 p , depending 
on the parameter < A < oo: 



±(v 2 h) dv, 



° x ( b ) = 6 E / f£e±(«i)<fo. 
± J 

V x {b) = bJ^J t>2»tQ±(vi)dv, 
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and a constant depending on A: 



where, for A = we define 



(v\) dv dx, 



(vi) dv dx. 



Finally, we derive formulas for the adjoint operators of B x , C x and V x . We 
begin with (B x ) • Let h,k G L 2 p . To simplify notation in this calculation, we 
drop the summation over ±. All calculations work similarly with the proper 
definition that includes the ±. 



(B x h,k) L2 = 



(j, p dv^j h(x) + j /j, e Q x (v 



k(x) dx 



fi p dv j h(x) + J fi 
Hp dv I h(x)k(x 



Xe As V 2 (s)h(X(s)) ds ) dv 
o 



k(x) dx 



\e Xs v 2 h{x)k{X{-s)) ds dv 



dx 



dv) k(x) + j fi e v 2 ( / Ae As fc( 



h(x) dx 



Therefore (with the ±) 



,k dv. 



(B X Y k = (W 4 dv\k + ^J ^v 2 Q> 
The computation of (C A ) is much simpler. Let 6 G K and k G Lp. One has 
{C x b,k) Ll =bj^ (j^jvtQi («i) dvj k{x) dx, 



and, therefore 



{C X )*k = J2j^ J f4Q±{vx)k{x)dvdx. 



The derivation of (T> j is similar and yields 
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2.4 The Matrix Operator M x 

In this section we rewrite Maxwell's Equations (|2.7[) in terms of the the un- 
knowns 4>, ip and b. The various operators acting on these unknowns are pre- 
cisely A X ,A X ,B X ,C X and l x . We show that Maxwell's Equations reduce to a 
simple, selfadjoint matrix operator A4 X , and that there exists a nontrivial solu- 
tion to Maxwell's Equations if and only if J\4 X has a nontrivial kernel. Note that 
in what follows, the dependence upon enters through the right-hand-side of 

(23). 

1. Rewriting Gauss' Equation (|2.7c[) as d 2 (j> + p = 0, we substitute (|2.12[) , to 

get 



= d^ + p = d^ + J(f + -r)dv 

= + I (^0(x)+/i^(^)-MeQ±^ + ^Q±(«2V') + ^Q±«l) dv 

± J 

= d^+lj2j^dvU(x)-J2j^Q±^dv 



(E J ^ dvU(x)+J2 J l*tQ±(W) dv + bJ2j 'ntQ 



±vi dv 



= -A x (t) + B x ijj + C x b. 

Thus, our first relation is 

- A x 4> + B x i> + C x b = 0. (2.17) 

2. Rewriting the second of Ampere's Equations (|2.7bj) as d 2 ip — \ ip + 32 = 
and repeating the same procedure, we have 

= d^-\^ + j 2 

= d 2 Jj-X 2 ^ + J2 I v 2 (txt()>(x) + ^(x)- fifQi^ + fitQ^ih^ + b^Qiv!) dv 
± •* 

= d 2 x 4> ~ A 2 y; + fe/ ^4 dv^ + Y,J v 2 »tQ±(v2i>) dv 

+ (^T J v 2 p± dv^j (j){x) Q±<i>dv + bJ2 J V2pf Q±vi dv 

= -A x i>- (B x )*(f> + V x b. 
where in the last equality we used (|2 . 1 3[) . Our second result is 

(B x )* 4> + A x 4' - V x b = 0. (2.18) 
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3. Finally we consider the first of Ampere's Equations (|2.7a|l : = XEi +jx = 
—\d x 4> — X 2 b + j\. Integrating over one period, we get 

1 f p 
X 2 b = — / jx dx. 

" Jo 

Plugging in the expression for the particle density /+ — /~ as before, we 
have 



\ 2 b 



pL hdx= p 



vi (/ + - / ) dv dx 



l4Q±(v 2 i>) + bfxfQ^vx) dv dx 



The first two terms vanish since fi is even in vx- Thus we are left with 

>?b=?pJ2j Q + dvdx. 

Denoting these three terms /, II, III respectively, wc find that 



P^J Q J v^tQl^dvdx 
bJ2 f I I ^e Xs v lf if(l)(X ± (s)) ds dv dx 
^E^/ f te^i-s^t^dsdvdx 
f [ [ Xe Xs V 1 ± {s)fif(l)(x)dsdvdx 

P , Jo J J-oo 



P ± 

= p(c A )*0, 

where for the third equality we changed variables (x, v) — > (X(—s), V(— s)) 
and used the facts that this change of variables has Jacobian =1 and that 
H is invariant under D, and for the fourth equality we changed the variable 
v — > — v and used the change of coordinates prescribed in (|2.16|) . Similarly, 

and, finally, 

III = bl x 
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by definition. Summarizing, we have 

{C x )* 4>-{V x )*^-Pb{X 2 -l x )=Q. 



(2.19) 



Motivated by the three equations (|2.17j) - (|2.19|) depending upon the parameter 
A, we define the matrix operator M x : Lp x Lp xl-> Lp x L 2 P x R 

/ -A x B x C x \ 

7W A = (B x )* A x -V x 

V \c x y -(v x y -p(x 2 -i x ) ) 

with domain Hp x ijj, x R. Formally, to prove our main theorem, it suffices 
to show that M. x has a nontrivial kernel for some < A < oo. Accordingly, we 
also define 

/ -A\ \ 
A4° = A° 2 . 
\ OP/ / 

Remark 2.6. As mentioned before, since <j> only matters up to a constant, 
we restrict the domain of M. x and M° to Hp x ifj> x R. Indeed, making 
this restriction is important. Letting (cf),ip,b) = uj riv = (1,0,0) we notice that 
M Utriv = for any A > 0. However, Utriv is a trivial solution that is of no 
interest for us, since it would generate a trivial solution (f,E,B) = (0,0,0). 
Moreover, multiples of utriv are the only trivial solutions. Indeed, whenever 
either -0 or b are nonzero, the linearized equations become nontrivial. 



3 Behavior for Small and Large A 

This section is the heart of this paper. To show that M. x has a nontrivial kernel 
for some < A < oo , we wish to show that there is some eigenvalue that crosses 
as A increases, starting at 0. To achieve this, we must somehow be able to 
keep track of the spectrum. Since we are working on the finite interval [0, P] we 
do not have to worry about a continuous spectrum. However, since A x and A x 
appear with opposite signs in A4 A , we expect to find eigenvalues near — oo as 
well as near +oo. Thus, to facilitate the counting, we truncate the spectrum, 
and leave only a finite number of eigenvalues. The following projection is defined 
so as to preserve the spectral properties of the two operators appearing along 
the diagonal: A® and A®- 

We postpone discussing the properties of the various aforementioned opera- 
tors until $71 However, these properties are used throughout this section. 

Consider the eigenvalues ot\ < oti < ■ ■ ■ of in ascending order, counting 
multiplicity. Let £j be unit eigenvectors associated to on chosen to be mutually 
orthonormal. We define the projection P n : Lp Q R™ to be: 

P n u={(u,& L A n , (3.1) 
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and, thus, P* : W l -> L 2 P is given by 

n 

P> = ^a^i, (3.2) 

i=l 

where a = (oi, aa, . . . , a n ) £ R™. Hence, 

iSP»u = (3-3) 

i=l 

is the projection onto the eigenspace associated with the first n eigenvalues of 
A° v ^ ^ 

Similarly, we define Q n to be the projection operator onto the eigenspace 
spanned by the n eigenvectors £1 of the operator A® associated to its first 

n eigenvalues /3i < /?2 < • • • < /?«. ■ We are now ready to define the approximate 
matrix operator M x : 

/ -P n A x P* P n B x Q* n P n C x \ 

M x n =\ Q n (B x )*P: Q n A x Q* n -Q n V x . (3.4) 
V (C X )*P: -{V x )*Ql -P(\*-l x ) J 

When A = 0, this definition reduces to 

/ -P n A\P* \ 

M° n = QnAlQ* n . (3.5) 

V o Pl° J 

Both matrices are finite-dimensional mappings R™ x E" x E — > E" x 1" x I. 
Let us list a few facts which will be useful for us later on: 

1. For all A > and for any n G N, .M^ is sclfadjoint. 

2. Let S(^) denote the spectrum of the operator A. Then S(M°) C S(M°). 

Proof. This is clearly true due to the properties of P„ and Q n , and the 
diagonal structure of both matrix operators. □ 

3. There exists some a* < such that [a*, 0) n S(7W°) = 0. 

Proof. Since the spectrum of A4 X is discrete (with no finite accumulation 
points) for all A > 0, there exists some oq < that is the greatest negative 
eigenvalue of M°. We can choose a* £ (co,0) freely. □ 



4. For fixed n > 0, Ai x varies continuously in A as a mapping '. 
E" x E" x R, for all A > 0. 
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Proof. For A > this is clear, since M x varies continuously by Lemma IT^l 
Thus we focus on the case A = 0. Since all norms on M. k are equivalent, 
it will suffice to check that M x — > M„ strongly as A — s- 0. Let u T — 
(4>,ip,b) e W x M" x M. We need to show that \\(M X - M°)u\\ ^ as 
A ->• 0. Since 



(M x -M° n )u = 



-P n {A x - A1)PZ 

Qn (B x yp* 

(C X )*P* 



PnB X Q* n 

Qn(A 2 — A 2 )Q n 
-(v x )*Q* n 



PnC X 
-QnV X 

-P(X 2 -l x + l°) 



it suffices to show that terms of the form \P n (A x -Aj)P*<f)\ or \P n B x Q* n ip\ 
tend to for fixed n, as A — > 0. This is clearly true due to the properties 
of the various operators discussed in §TJ □ 

5. For fixed n > 0, PnPn is bounded both in Hp and in Hp. 

Proof. This follows from the fact that P^Pn maps Lp into Hp and the 
closed graph theorem. □ 



Lemma 3.1. For any g £ Hp, P*P n g 
Similarly, Q* n Q n g -> g in Hj, and in H 2 P . 



g in Hp and in Hp, as n — > oo. 



Proof. We prove for P^Pn', the proof for Q* n Q n is similar. Assume that g = 
Y^jLi9j£j m ipi where £j are eigenvectors of A\, as defined above. Since 
g € i?f> = D(^4j), we know that A\g € £p and, therefore, there exist some 
/Hj such that A\g = Y^jLifij^j w i tn Yl'jLi \Pj\ 2 < 00 • ^ n f act ; by taking the 
Lp-inner product of A^g with we easily see that = gkO-k, where is the 
fcth eigenvalue of A®. 

Using Poincare's Inequality (twice) and then the triangle inequality we have 



OO 




OO 




OO 




OO 




< c 




< c 




+ c 




j=n 






^ 2 P 


3=n 







The first term tends to since A\ + d 2 is a bounded operator on Lp (see 
Lemma l7.2t |T))') . and, therefore, this term is controlled by C Y^jL n 9j£ 

oo since g £ Lp. The second term tends to since 



which tends to as n 

^ <} y.j ,. 9j€j = „ and, as mentioned above, Y^Li \9j<*j\ < °°. 

This shows strong convergence in Hp. The strong convergence in Hp is due to 
interpolation between Lp and Hp. □ 



Corollary 3.2. Considering the restrictions of P^P n to Hp, the lemma implies 
by the uniform boundedncss theorem that 



sup || P*P n || 



(3.6) 
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Lemma 3.3. The operator P*P n can be extended to H p . The sequence of 
extended operators {PnPn}^ = \ converges strongly to the identity. 

Proof. Our basic tool in extending the domain of definition is the canonical 
identification of Lp with a subspace of Hp 1 via the scalar product. With that 
standard definition at hand, we can now extend P„P n to Hp . For brevity, 
denote J„ = P*P n '■ Lp — > Hp. As noted above, the restriction to J„ : Hp — > 
Hp is bounded. Hence the dual operator J* : Hp 1 — > Hp 1 is bounded with 
the same operator norm (see [26) III-§3.3] ) . In fact, J* is an extension of J„ to 
Hp . Indeed, if g £ Lp and / 6 Hp, and using the (Hp, Hp 1 } paring, one has 

(/, J* n g) = (Jnf, g) = (J n f, g) L % = (/> J ng) L % • ( 3 - 7 ) 

Finally, we must show that the sequence {Jn}^=i converges strongly to the 

identity in Hp 1 as n — > oo. In view of (|3.6p . also sup n 1 1 «^ T * 1 1 ^ 1 > ^ ^ 1 < 00 • 

Therefore, it suffices to prove that J*g — > g for g in a dense subset of f/p" 1 . Since 
we already know that strong convergence holds for g 6 Lp, wc are done. □ 

Lemma 3.4. Let = (c/) n ,ip n ,b n ) <G R"xR™xR, and suppose that (P*4> n ,Qn' l Pn,b n ) — > 
u a = (<^OiV>Oi&o) strongly in Lp x Lp x R as n — > oo with P*^ n ,Q^tjj n uni- 
formly bounded in Hp. Suppose that A n — > Ao € [0,oo) as n — > oo. Then 
A4 Xn u n —> Ai x °uo weakly, in the sense that 

( M nUn,V n ) RnxWnxM ^(M X °U ,v) L2pXL%xM . 

for any v G L P xLpXl with v n € R™ x R™ x R being the projections of the 
two first corrdinatcs of v and the identity on the third. 

Proof. Let us write Ai^u^ and Ai x °uo explicitly: 



PnA 1 " P„ 
M^Un = | Q n (B X ")* P* 

(c^yp* 

-A Xo B x « 
M x "u = I (B x °y A x " 

\c x «y {v x °y 

We show the weak convergence term by term, by its location in the matrix. 
Let us demonstrate for the first term - the rest are either similar, or simpler. 

We want to show that P n A Xn P*4> n —> A l °4>° weakly. Since the sequence 
A Xn P* 4> n is uniformly bounded in L p by assumption, it suffices to test the weak 
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continuity by taking g £ Hp which is dense in Lp. We have 



P n Al"P*<i> n ,P n g) -{Al <t>°,9 



< 



P*P n A^P:4> n -Al°<tP,g 



P* n P n [A^ -A\°)P*<t> n ,9 



{P*P n -I) A^P^g 
(je?{P* n ct> n -ctP),g) L2 



I + 11 + III. 



Terms / and II are treated very similarly: 



p: 



■^1™ -^1° ) PnPn9 



A„ 



<ll^ 



Ay" Ai° ) P n P n g 



and 



11 = 



P^ n ,At° {P*P n -I)g 



< \\P, 



A? {P*P n -I)g 



->0, 



since H-P^raHj^ is bounded, and since A\ n - A]° tends strongly to as an 
operator Hp — > Lp and PnP n tends strongly to / in Hp by Lemma |3~T1 Finally 



III 



p: 



iO \ An 



,A?g 



since P* 



b° strongly in Lp by assumption. 



□ 



Lemma 3.5. For n sufficiently large, has exactly K n := n — neg{A\) + 
neg(A^) + neg(l°) negative eigenvalues. 

Proof. Since A4^ is diagonal, we may consider each entry along the diagonal 
separately. The number of negative eigenvalues of — P n A\P^ equals the number 
of positive eigenvalues of P n AiP*, namely pos{P n A\P*). Since we assume that 
kcr.A" = {constants}, and since our domain does not include constant functions, 
P n AiP* has no null space. Therefore, pos{P n A\P* L ) — n-neg(P n AiP*). How- 
ever, for n that is sufficiently large neg(P n A\P* l ) = neg(A\) by our definition of 
P n . Thus, the contribution from the first term on the diagonal is n — neg{A\). 
The contribution to the negative spectrum from the next term on the diagonal is 
n&g{J\^) = negiQnA^Qn) f° r n that is sufficiently large. The contribution from 
the last term depends upon the sign of I and is denoted by neg(l°). Combining 
these observations we get our result. □ 

Proposition 3.6. There exist N £ N and A* > such that for any n > N and 

for all A G [0, A*], neg{M*) > neg(M° n ) = K n . 
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Proof. We show the existence of A* and N by contradiction: If they do not exist, 
then for any k > 1 there exist n k > k and \' k < £ for which neg(A4 n k k ) < K nk . 
We choose fc large enough for Lemma [3.51 to hold, and fix it. Since M x k is a 
continuous mapping in A, its spectrum varies continuously with A. Therefore, 
since neg{M.^ lk ) = K nk > neg{Mn k k ), at least one eigenvalue in the negative 
part of T,(M X ) must cross from left to right as A varies from to X' k . Since 
there is no spectrum of M^ k on the interval (er*,0), this eigenvalue must also 
cross this interval, and, in particular, it must cross a* at some value A& € (0, X' k ). 

To summarize this argument, our contradiction asserts that there exist n k — > 
oo and Afc — > for which 

M%u nk =a*u„ k , (3.10) 

where ^ uj lk = (</>„ fc , tp nk , b Uk ) £ 1" x I" x i To simplify notation we drop 
the "fc" index and simply write A„ instead of A& and n instead of n k - Our plan 
is to show that this contradicts the fact that a* is not an eigenvalue of A4°, 
as follows: First we show that u n has some nontrivial limit in an appropriate 
space; then we use Lemma I3~4l to show that the operators M Xn converge (in 
the weak sense defined in the statement of that lemma) to . 

Since P* and Q* n both have trivial kernels, we may normalize the eigenvectors 
as follows 

\\PMl% + \\QM\l% + K\ = 1, (3.11) 

and, therefore, there exist the two weak limits in L 2 P and the limit in R, respec- 
tively (after extracting a subsequence) 

P*<f> n ^<f>, Q^ip n ^ip, b n ^b. (3-12) 

Our goal is to show that {<f>, tp, b) is nontrivial, and that 



(3.13) 



thus reaching a contradiction to the fact that a* ^ ker(A^°). Note that we do 
not have to worry about showing that (0, ip, b) is not a multiple of (1, 0, 0), since 
(1, 0, 0) is in the kernel of M. x for all A > 0, but a* =/= 0. We begin by showing 
that both P*4> n and P^ipn are bounded in Hp. The first row of (|3.10j) is 





P 



(-A\ n Pl4> n + B x -Q*J> n + C x -b^ = a*4>n- (3.14) 



Write 4> n = (</> r \, 0^, ... , </>"), and take the inner product in 1" of (|3.14[) with 

4>n- 

0-*\4>n\ 2 = 0„-P„(-^"P„*0„+S A "Q;Vn+C A "6„) 

= (P* K , -A\~ P* <t>n) L2 + K ■ Pn {B X - Q*jM + C A " b n ) =1 + 11. 
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Since 1 > \\P*(j)n\\ 2 L 2 p = £fe=i |C| 2 = |0„| 2 , the left hand side of the above 
equation is uniformly (in n) bounded. Moreover, 



= (P^dlP^n), 



-\\d x P* n K\\ Ll +h+h 
where it is easily seen that 



I'll 



p: 



E / Me fo\ K 



< IIP, 



E SU P / iMe | cfo < C < oo 



and, as in the proof of Lemma l2"31f Tj) . we have 



\h\ = 



£ P*<j> n E / Me Q± (O^ 



E 



\ An* 



P*<t>n{x) P*4>n{X ± { S ))dvdxds 



< 2 ||P! 



£<C||OX <<7<oo. 



In addition, 
l/JI 



P„ U?" 



(P:<j )n ,B x -Q* n yj n +C x "b n ) l 



3 



\Qn"<P- 



< C < oo. 

by Lemma 17.21 Thus ||9 2; P*^i n || i 2 is bounded uniformly in n in Lp, so that 
||P^^n|lij| is uniformly bounded. Therefore P*<j) n converges strongly in Lip and 
weakly in Hp to <f> by (|3.12[) . Similarly, we show that Q* m %p n is bounded in Hp 
by considering the second row of (|3.10[) : 

(3.15) 



Q m (S A ") P*<j> n + A^P^ n - V x ~b n ) = a*ip n 



The analysis is similar. Hence Pn4>n and PjtV'n have strong limits in Lp and 
weak limits in Hp , which must be <j> and ip respectively. Due to p. lip the limit 
{P*4> n ,P*ip ni b n ) cannot be trivial: {^,tp,b) ^ (0,0,0). 

It remains to be shown that, in fact, <p and tp lie in Hp and hence, in the 
domain of A4°. Consider again equation l3.14l We rewrite this equation, keeping 
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only the Laplacian on the left-hand-side, moving all other terms to the right- 
hand-side and applying P*. We get 

KPndl (P„>„) = - P^n IT, j £ P n<t>n + P*P n £ / Me Q± (KM dv 

- P:P n B x -Q* n ^ n - P*P n C x «b n + a*P*j> n . 

The last term may be written as a* P*P n P^(t> n , so that we may denote the 
entire right-hand-side as P*P n h n for brevity. Now, we know that h n £ L 2 p has 
a strong limit in Lip, and P*P n —tI strongly in Lip. Therefore, the right-hand- 
side converges in Lip. Using Lemma |3~31 and since P^n converges weakly in Hp 
to cf>, the left hand side converges weakly in the Hp 1 sense to d%</>. By elliptic 
regularity one can bootstrap and deduce that, in fact, <f> £ Hp. 

Finally, using Lemma 13.41 we know that the approximate equations (|3.10p 
tend (weakly, in the sense of the lemma) to the equations M°u = <j*u. But 
a* £ E(M°). This contradiction ends the proof. □ 

Proposition 3.7. If the null space of A® is trivial, then for any n > N and 

for all A £ [0,A*], neg(A4 x .) = neg(A4„) = K n , where N and A* are as in 
Proposition 13.61 

Proof. Using the same contradiction argument as in Proposition 13.61 we can 
show that pos(M.n) > pos(M^) for small A and large n. We conclude that an 
increase in either the number of negative eigenvalues or positive eigenvalues can 
only be due to zero eigenvalues of M„ that move left or right as A increases. 
However, under the assumption that the null space of A% is trivial, M? n has a 
trivial kernel as well, and, therefore both neg(A4 x ) and pos(A4 x ) must remain 
constant (and equal to their values when A = 0) for small A and large n. □ 

Lemma 3.8. There exists A* > A (A as in Lemma I7.1[) such that for every 
n £ N and for any A > A*, M x has precisely n + 1 negative eigenvalues. 

Proof. Since M x : R n x R" x R ->• R n x R" x R is symmetric, it has 2n + 1 
eigenvalues, all real. Letting ip £ R™, we have 




M X n ( J& j , ) ) = (QnA^Qn^^)^ = (A Q*J>,Q*J>) L > > 

(3.16) 

for all A > A by Lemma 17.11 This implies that Ai x is positive definite on a 
subspacc of dimension n, and, therefore it has at least n positive eigenvalues. 
Similarly, we now show that there exists a subspace of dimension n + 1 on which 
M x is negative definite: Let (</>, 0, b) £ R" x R" x R and consider 




= - <^P>, Ptf) L% + 2 {c\ P*4>) Ll P(X 2 i> 

(3.17) 
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We estimate the middle term as follows: 



(c x b,p: 



\C x b\ 



<2 C x b L2 \\P* 



< 



e 2 \\p: 



Letting e 2 = i, we have 




-P(A^-Z A )fe 2 + A||C A & 



A 



Using the fact that ,A A > 7 > 0, for all A > A (see Lemma [7T][6])), this 
expression is negative for all <j> G R™ and b G R, since Z A and C x are both 
uniformly bounded. Therefore, there exists a A* > such that for every A > A* 
there exists an n + 1 dimensional subspace on which M x is negative definite. 
We conclude that 



neg (Ml 



1. 



for all A > A* 



Notice that A* does not depend upon n. 



(3.18) 
□ 



4 Limit as n — > oo 

Lemma 4.1. Let A*,A*,A^ be as above. Fix any n > N. Then there exists 
A„ G [A*, A*] such that M Xn has a nontrivial kernel. 

Proof. As we have seen above, A* and A* do not depend on n. We apply a 
simple continuity argument: M x is continuous in A for each (fixed) n in the 
sense that if a > 0, then there exist C, 6 > such that 

l|M A -X ri J| <C\X-a\ 

for A G (0,oo) and | A — | < 5. This follows from Lemma \7. 41 By Proposition 
13.61 M x has at least n — neg (Af) + neg (A®) + neg (Z°) negative eigenvalues. 
By Lemma 13.81 M^ has exactly n + 1 negative eigenvalues. Since M x is a 
finite-dimensional operator, its set of eigenvalues varies continuously with A. 
Thus, if 

n - neg (A") + neg (A%) + neg (Z°) > n + 1 (4.1) 

then at least one eigenvalue must cross for some A n G (A*, A*). In particular, 
the A value for which is a (nontrivial) eigenvalue, has a corresponding (non- 
trivial) eigenspace. Unraveling condition (|4.ip , we get the equivalent criterion: 

neg (A° 2 ) > neg (A") + neg(-l°) (4.2) 

which is precisely the main assumption of Thcorcm[TJ In the context of Theorem 
[5] one would invoke Proposition 13.71 instead of invoking Proposition ^. 61 and get 
the criterion 

neg (A° 2 ) ? neg (A ,) + neg(-l°). (4.3) 

□ 
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Lemma 4.2. There exists < Xq < oo and a nontrivial Uq = (4>o,ipo, bo) that 
is not a multiple of (1,0,0) such that 

M Xo u o = 0. (4.4) 

Proof. This proof is very similar to the proof of Proposition ^. 61 We have seen 
that for all n > N, A4 X has a nontrivial kernel when A = A n if (|4.1[) is satisfied. 
We show that (|4.4[) is satisfied weakly, with uq extracted by some compactness 
argument from nontrivial elements in the kernel of VW^™, denoted by u n . As we 
have seen in Lemma 14.11 the equation 



/ -PnA^P* PnB X "Q* n P n C X " \ ( (f> n \ / 

M x "u n =\ Q n (B x ")*P* Q n A Xn Ql -Q n V x - = 

V (c x »y P * q* -p(\i-i x ~) )\b n ) \ o 

(4.5) 

has a nontrivial solution for all n > N. Here < A* < A„ < A* < oo. Let us 
extract a subsequence A„ — > Aq. Wc want to show that the "limiting" equation 



/ -A Xo B x « C x « \ ( cj> \ ( \ 

M x "u = (B X «Y A Xa -V x « ^o = (4.6) 

V (c x »y -(v x «y -p(xi-i x ») J \ b J \o J 

is satisfied nontrivially. We follow the procedure of Proposition 13.61 showing 
that {P*<p n ,Q* n ^n,b n ) = ul -> t# = (0o, Vd. M + (0, 0, 0) in ffj, x J, x K: 

1. Wc normalize the vectors u n as in (|3. 1 1|) : 

\\P*<t>n\\Ll + \\QML% + \bn\ = l- (4.7) 

2. We take the inner product of the first row of (|4.5p with <j> n , to obtain the 
equation 

K ■ P n (-A\ n P*<t> n + B Xn Q* n ip n + C A "6„) = 0. (4.8) 

Showing uniform boundedness of P* x <i> n in Hp is identical to the calcula- 
tions performed in the lines following (|3.14j) : We show that the Lip norm 
of all terms in (|4.8|) is uniformly bounded in n (except for the Laplacian), 
and, therefore, by integrating by parts we obtain the uniform Hp bound. 
We conclude that P^</> n converges in Lip. 

3. Using the second row of (|4.5[) we show that Q* n ii> n converges in Lip. 

4. We bootstrap our convergence problem, showing that, in fact, P*4> n con- 
verges to (f> in Hp, by applying P* to the first row of (|4.5[) . Similarly, to 
show that Q* tp n — s- tp in Hp we apply Q* to the second row of (|4.5[) . 
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Moreover, since P^n G Hp all have mean 0, so does fa. Therefore u is not a 
multiple of (1,0,0). The fact that the terms in the equations (|4.5[) tend to the 
terms in the equations (|4.6[) weakly follows from Lemma \'3. 41 This finishes the 
proof. 

□ 

5 Construction of a Growing Mode 

We finish the proof of Theorem [T] by verifying that the nontrivial element Uq 
that we found above satisfies the linearized RVM System. For ease of notation, 
we drop the "0" subscript, so that we simply have u T = (cf>,ip,b), and A. The 
equation we verified in the previous section is 

( -M & c x \ / cb \ 

M x u = (B x )* A\ -V x V =0, (5.1) 

V (C X Y -(V X Y -P(X 2 -l x ) I \ b j 

with u nontrivial and not a multiple of (1,0,0), and where < A < oo. We 
begin by defining f ± (x 1 v): 

f ± (x,v)=± l it<t>(x)±4iP(x)T»i [Qi^-Qii^-bQlv!] . (5.2) 
In addition, we define 

Ei = -d x <j> -Xb E 2 = -\ip B = d x ip 

and 

P = J(.f + - D dv ji = J Vi(f+ - f~ ) dv, i = 1, 2. 

Lemma 5.1. Gauss' Equation (|2.7c[) holds. 
Proof. We use the first row of (|5.1I) . 

d x Ei = -d 2 x ct> 

= dv(j>- J ntQ±<j>dv + Jnpdvip + J fitQ±(v 2 ip)dv + J '^C±(«i) dv 6 J 

{f + -Ddv = p. 

□ 

Lemma 5.2. The Linearized Vlasov Equation (|2.5p holds. 
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Proof. We recall that the Linearized Vlasov Equation (|2.5[) is 

(d t + D*) /± = TufviEt ± »± Vl B T (Me «a + ) ^2. 

We let g G C* ([0, P] x M 2 ) be any test function. We show it for the electrons, 
/~, and drop all "— " superscripts. Showing that the linearized Vlasov Equation 
holds for /+ is identical. We write: 

f (Dg)f dvdx = [ [ (Dg) {-^{x) - H p ip{x) + \x e [Q x </> - Q x (v 2 yj) - bQ x v 1 ] ) dv dx 
o J Jo J 

= I + II + III + IV + V. 

For the terms I and II we use the fact that D is skew-adjoint and that 
fi is invariant under D, to "move" the operator D over from g to <fi and ip 
respectively. Thus we focus on the terms III, IV, V, where the definition of Q x 
is important, and the fact that (x,v) — > (X, V) has Jacobian = 1. 



III = J {Dg)ii e Q X (j)dv dx 



o 



,P , 

\e Xs / / fi E (Dg)(x,v)<f>(X(s; x,v)) dv dx ds 
-ao Jo J 
o r P r 

Xe Xs / / n e (Dg)(X(-s),V(-s))(j)(x) dv dx ds 

-oo JO J 

p r r° ( d \ 

H e / Xe Xs l-—g(X(-s),V(-s))]ds <f>{x) dv dx 

( i-O 

2 As 



[i e \ -Xg(x, v) + I X 2 e Xs g{X(-s),V(-s)) ds j tj>{x) dv dx 
pL e X(j){x) + He I X 2 e Xs (j)(X(s), V(s)) ds > g(x, v) dv dx 

J — oo J 

J J {-He4> + HeQ x 4 1 } g dv dx. 

IV = -A J J {-fi e V2ip + HeQ X {v 2 ^)} g dv dx 
J J {— b/i e vi + bfi e Q x vi} g dv dx. 



= X 
Similarly 



and 



V = —X 
Thus 
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r-P r pP 



(Dg) fdvdx = J J {p e D(j) + p p Dip} g dv dx 

+ A I [ {-^(1) + ^eQ X (t) + ^eV2lp - fJ-eQ X (v2^) + bfl e V! - b^i e Q x vi} g dv dx 

A {-Me</> - MpV' + Me [Q A - Q A {V21P) - 6Q A «i] }gdvdx 

{XfXptp + n e D<j> + /i p Dip + Xfx e V2tp + Xbp e vi} g dv dx 
{A (/ + [ipip) + fJ, e D(f> + i-LpD^j + Xp e v 2 ip + Xbii e vi} g dv dx 



Therefore, weakly, / satisfies the equation 

(X + D) f = — fi e Dcj) — pipDi/j — Xfipip — Xp e v 2 %jj — Xbp e v\ 

= -fj, e vid x <j) - (ipVidxi/j + fi p E 2 - Xp e v 2 Tp - Xb^ e vi 
= HeV\E x - [ipV^B + (pp + p e v 2 ) E 2l 

which is precisely (|2.5p . □ 
Lemma 5.3 (Continuity Equation). The relation d x ji + Xp = holds. 
Proof. Integrating the Linearized Vlasov Equation with respect to v, we get 

J (d t + D*) f ± dv = J (T^VxEi ± ViB t (^h + /x±) E 2 ) dv = 0, 

where we use the facts that /i* is even in vi, and that d\i^ jdv 2 = \i^rv 2 + /i^ 
is a perfect derivative. Subtracting these two equations, replacing the time 
derivative by a factor of A, and using the fact that D consist of three terms, 
of which only the first is not a Vi derivative, we have: 

= f ( A + D± ) / ± dv = ^ + 9^1 vif* dv = \p + d xJ1 . 

□ 

Lemma 5.4. Ampere's Equations (|2.7aj) and (|2.7b[) hold. 

Proof. We first want to show that 

XEi = -h. (5.3) 
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Recalling the definition E\ = —d x <\> — Xb, we wish to show that 

\ 2 b=-\d x cf> + j 1 . (5.4) 



Let us show equality of the derivative with respect to x, and then equality 
of the integral with respect to x. Differentiating once, and using the continuity 
equation, we get 



= -dl<P + -d x]l 
= -d 2 A - p 

± J 

= A X (f>-B x ijj-C x b 

which is precisely the first row of (|5.f [) . This verifies that the derivatives are 
the same. 

Now we turn to the integral of (|5.3[) : Plugging in the the relation A 2 6 = 
js (C x )* 4>-j; (£> A )* ip + l X b (which is obtained from the last row of (|5.ip) into 
(|5.4p . it suffices to show 

-\d x <f> + 3i = ^ (C A ) *4>~j (V X ) * ^ + l X b- 

Writing in detail the expressions for j\ and for the operators on the right hand 
side, we need to show that 

-\d x <f> + [vi {^<f> + 4 $ - ^ [Q±<i> - G± (flaVO - 6C±«i]) dv 

= T^l J [^2i(«i)'/ , -«2M e ± Q±(«i)V' + 6«iMc ± Qi(«i)] dvdx. 

Since y, is even in v\ , we may drop the first two terms in the integral on the left 
hand side. Therefore we need to show that 

-\d x <t> = I »ia£ [G±^ ~ G± (^VO - b Q± dv 

+ I ^ [Q±(^i)^-^Qi(vi)'ip + bv 1 Q x z (v 1 )}dvdx. 
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Integrating this equation over the period P we get 



= JlJ o J <W [Q±<t> - Q± («a^) - («i)] <*« 

+ / t4[Q±{^)4>-V2Qi{vx)^ + bvxQl{v x )\dvdx 

, Jo J j- Jo J 



which indeed holds due to the change of variables (x,v) — > (X,V). Therefore, 
both the derivatives and the integrals are equal. Hence (|2.7a[) holds. 

Now we turn to show that the equation XE2 + d x B = —32 holds. We again 
recall our definitions 



E 2 = -\^, B = d x i/>, 32 = J v 2 (/+-/") do, 
that imply that we need to show 

-aV + ^ = - J Mf + - r) dv 

But this is precisely the second row of (|5.1|) . 

□ 

Our proof of Theorem [T] is now complete. 

Proof of Theorem [H The proof of Theorem [2] is identical to the proof of Theo- 
rem [TJ with the only difference being that in Lemma |4. II we invoke Proposition 
13.71 instead of invoking Proposition 13.61 □ 

6 Examples 

6.1 Homogeneous Example 

We start with the simple homogeneous case, where there is no x dependence. 
This case is so simple, that certain properties of the operators can be calculated 
explicitly. In this case e ± = (v) ,p ± = V2- Hence = fj,~ , and we can therefore 
drop the ± in this example. 



Lemma 6.1. In the homogeneous case, 



1 r P 



V ["i(v)h{x)] = j(v) ■ — I h(x) dx. 
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Proof. This is straightforward, since D reduces to the simple differential opera- 
tor v\d x in the homogeneous case, since there is equilibrium magnetic field. □ 

Proposition 6.2. 1. Any homogeneous equilibrium that satisfies 

(a) / (i e iif dv < 0. 

(b) / /z e dv < 0. 

(c) J \i c v\ dv — J fi p V2 dv < 0. 

is unstable. 

2. There exists such an equilibrium. 

Proof. 1. We verify that the above conditions verify the conditions of the 
first part of Theorem [T] which would imply that the solution is unstable. 
Our plan is to show that the conditions listed in the proposition guarantee 
1° < 0,neg (A%) > 1, and neg (Al) = 0. Since \i has no x dependence 

1° = — J J [i e v\V (v\) dv dx = J [i e v\V [%\) dv = J \i c v\ dv < 0. 
Next, we have A\h — —d%h — (/ /i e dv) h+ -p J fi e J Q P h dx dv, so that 
(A%h) = {h'f dx - (J fi e dv^j 




> 



by Holder's Inequality since J fj, e dv < 0. Thus neg(Ai) = 0- As for A°, 
we write 

A\h = —d^.h — V2[i p dv^j h + — J v\\i e dv J h dx. 
Therefore, choosing h = 1, we have 

(All, 1) = -P J v 2 Hp dv + P J vlfi e dv < 0, 
so that neg(A 2 y ) > 1. 

This verifies that the conditions imply that neg{A2) > 1 > = neg(A\), 
and that 1° < 0. Therefore, by Theorem [T] such an equilibrium is unstable. 

2. We construct an explicit example. Let /i(e,p) = a{e) with a(e) = 7(e) + 
77(e), where 7(e) = e — 1 on the interval [1, 2) and otherwise, and 77(e) = 
exp [— (e — 2) 2 ] on [2, 00) and otherwise. Even though a(e) is not smooth, 
it can be approximated by a smooth function that will still verify the 
calculations below. We verify the conditions of the proposition one by 
one: 
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(a) We need to verify that J [i e v\ dv < 0: 

\j, e v\ dv = a' v\ dv — / (7' + rf) v\ dv. 



We make the change of variables V1V2 — > r0, so that r 2 = vf + v\ 
e 2 — 1, to get 



jJL e V\ dv 



(V +rf)—i ^-rdrdd 

1 + r z 



Jo ' 



cos^ de 



1 + H 



73 1 + 



dr 



Now: 



rv^ 3 

J = / ? dr ps 0.6 

1 + r^ 



dr 



which, indeed, verifies the first condition 
(b) We verify that J fj, e dv < 0. 



2 (y/l + - 2) cxp - (Vl + »- 2 - 2^' 



1 + r 2 



dr ps -2.5 



/i e di> = 



a' dv 

271 



d<9 





/>oo 


/ r dr + 


/ rf r dr 


Jo 


JV3 



= 2tt 
ps 2tt 



o /"oo r . . 

--y 2 (Vl + r 2 - 2J cxp - ( sjl + r 2 - 2 J 



'V3 



dr 



< 0. 



(c) We have J /Lt e u| dv — J [i v v>i dv = J ii e v\ dv since does not depend 
on p. The proof that this integral is negative is identical to the above 
proof that J /J, e i 2 dv < 0. 



□ 



Thus our assumptions are all verified. This implies instability. 
6.2 Weak Magnetic Field 

Proposition 6.3. There exists an inhomogeneous, nonmonotone, purely mag- 
netic equilibrium that is unstable. 
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Our goal is to construct an explicit purely magnetic equilibrium, for which 
we can conclude, using our main result, that it is unstable. Again, we would 
like to verify the conditions of Theorem [TJ 

1° < and neg (J%) > neg (.4?) . 

Therefore, we construct an equilibrium /i ± (e,p) for which 1° < 0, neg (A®) > 
1 and neg (A®) — 0. The main idea of the construction is to consider an 
equilibrium that is almost monotone. We separate x, v space into the sets Sg 
(good) and Sf 1 (bad) where are negative or positive respectively. Then we 
show that if S b are not too big (in measure) we are essentially in the previously- 
known monotone situation, and we can easily investigate the properties of the 
operators A\ and A%. To choose a magnetic potential we consider the ODE 
(|2.4|) for the magnetic potential, which can be written as 

d^° = 2 J 'v2ti-((v) ,v 2 -ip°(x)) dv (6.1) 

using a simple change of variables, similar to the ones demonstrated in the proof 
of Lemma I7.2t [2j). and recalling that /i + (e,p) = fj,~(e,—p). For simplicity, we 
write (|6.ip as 

In |34j it is shown that a sufficient condition for the existence of a purely 
magnetic equilibrium is the existence of magnetic potential that solves this ODE. 
Indeed, there exist periodic solutions. Moreover, assuming that u~(e,p) is even 
in p we readily see that g(0) = 0, and assuming that J v 2 (J-~ ((v) , v 2 ) dv > wc 
see that <7'(0) < 0. This implies that the origin is a center and that there exists 
eo > for which there is a family of periodic solutions ipe( x ) with < e < e 
and periods T^o depending upon e, satisfying (|6.ip . with 

1. \ip° I j -> as e ->• 0, 

2. T^o — > P cr as e — > 0. Here P cr is the period associated with the homoge- 
neous case tfj° = 0. 

By readjusting the starting point, wc may assume that obtains its minimum 
at and at T^o, that it has a single maximum, obtained at hT^o, that it is 
strictly increasing on (0, \T^o) and finally that it is symmetric with respect to 
^T^o. We now show that for e that is small enough, and with an appropriate 
choice of particle distribution fx~(e,p), this distribution is linearly unstable. As 
discussed above, wc arc flexible in our choice of particle distribution u~(e,p), 
as long as the following conditions are satisfied: 

1. /i~(e,p) is even in p 

2. / v 2 (i- ((v) ,v 2 ) dv>0 
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3. sup^ e and \Sb\ satisfy the relation 

SUP ^" < Sliwy (6 - 2) 

where |5f,| is the measure of the set Sb = > 0}. 

The first and third conditions are clearly easily satisfied with the right choice 
of fi~{e,p). As for the second condition: Since p is essentially V2 (they differ 
by ij>® which is a perturbation of 0), the second condition essentially states that 
/ pfjL~ dv > 0. This is satisfied since fi p is an odd function of p. 
We now show that A® is a nonnegative operator: Recall that 

A\h = -d 2 x h - (^J dv^jh + J ^Vh dv 
where V is the projection operator onto kerD. However, 

d = Vl d x - i> 2 B°d vl + wiB°a„ 2 , 

where B® = d x ip®. Thus V and A\ both depend on e. We denote them by 

€ 2 

V € and A{ . Let K e = 7^3- be the constant given by Poincare's Inequality, 

which depends on TLo and thus on e. This constant satisfies J h 2 dx < 

K~ x J (d x h) 2 dx. Since T^o depends on e continuously, so does K e . We let 

2 

Kq = -pr be the Poincare constant associated with the homogeneous case e = 0. 
We have the following: 

Claim. If (|6.2p holds then there exists some e' > such that A { e is nonnegative 
for all < e < e' . 



Proof. We first show for e = 0, and then conclude for small values of e by a 
certain continuity argument. For brevity wc drop the e. Letting u £ Hpo not 
identically 0, we have: 

(A\u,u) j 2 = / (d x u) 2 dx — \ // fi^ dv u 2 dx + // ^fV ± (u)udvdx 
p Jo ±JJsf ±JJsf 

-^2 ln±dvu 2 dx + Y^ I fitV ± (u) u dv dx 



± 

I - II + III - IV + V. 



We need to show that this is positive for all u. Since V ± are projection operators, 
we easily have \III\ < \II\ and \V\ < \IV\. Since II has a 'good' sign, we have 
that —II + III > 0. Terms IV and V have the 'wrong' sign. Using Poincare's 
inequality we can estimate: 



// ± £ dv u2 dx - H max/i^ )\ s b\ ' 
± •'•' s b ± \ s >> ) •'° 



u 2 dx < ^-K^ 1 I (d x u) 2 dx = i/. 
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Thus (A\u,u) >I-IV + V>\I>^ \\u\\ 2 L 2 p > 0. Finally, since K e varies 
continuously with e, we conclude that there must exist an e' as in the claim. □ 

Next, we turn to the operator A®- Again, one should denote A 2 ' e to make 
clear the dependence upon e. It is shown in section 4.3 of [34] that for e suf- 
ficiently small A% c has a negative eigenvalue if A^° has one. The same proof 
still holds in our situation. 

Our last duty is to show that = J 6 J v\[i~V t («i) dv dx is negative. 

By Lemma \2. 41 V preserves parity with respect to v\. Thus v\V € (vi) > 0. But 
since we chose fi~(e,p) to be such that the domain where fi~ > is negligible 
compared to the domain where [i~ < (and fi~ is bounded by a small positive 
number from above), we can clearly pick n~(e,p) such that 1° < 0. 

7 The Operators 

For the sake of completeness of this paper, we prove the important properties 
of our operators in full detail. These proofs appear in very similar form in 
[341 133J . One significant difference is that we cannot use fi e as a weight as it 
may vanish, and we therefore use the weight w introduced in the introduction. 
Another notable novel part is Lemma l7.1l[6|) concerning the positivity of A\ for 
large values of A, which has a rather lengthy proof. Most of these lemmas are 
consequences of the properties of Q± discussed in Lemma \2. 5 1 

Lemma 7.1 (Properties of Ai,A 2 )- Let < A < 00. 

1. Ai is sclfadjoint on L 2 PQ . A\ is selfadjoint on L 2 P . Their domains are 
Hp and Hp, respectively, and their spectra are discrete. 

2. For all h(x) <G Hp , \\A$h - A^h\\ L 2 — s- as A — s- 0. The same is true for 
A\ with h(x) G Hp. 

3. For i = 1, 2 and a > 0, it holds that \\Af - Af || = 0{\\ - a\) as A -> a, 
where |j • || is the operator norm from Hp to Lp in the case i = 1, and 
from Hp to Lp in the case i = 2. 

4. For all h(x) G Hp Q , \\A$h + d%h\\ L 2 p ->• as A -> 00. 

5. When thought of as acting on Hp (rather than Hp ), the null spaces of 
Ai and A\ both contain the constant functions. 

6. There exists 7 > such that there exists A > such that for all A > A, 
A\ > 7 > 0, i= 1,2. 

Proof. 1. Note that as mentioned in Remark 12.61 we are only interested in 
the action of A\,Ai on Lp and not Lp, but that does not matter for 
the purpose of this lemma. We first show that the perturbations of the 
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Laplacian in A$ are bounded operators for i = 1, 2 and all A > 0. For the 
case A > 0, a typical such perturbation may be estimated as follows: 



J £Q±hdv 




< (sup J \fj,f\dv 

< (sup J \i4\dvj \\Q±h 
<Y t c± \\^ h L^ c W h h% 



Q±h\ dv\dx 



1/2 




Q x ± h\ dv\dx 



1/2 



due to the decay assumption ()1.5[) and to the estimate ()2.15j) . Here, C is 
a universal constant, depending only on the equilibrium. The case A = 
is even simpler, since the operators Q± are replaced by the projection 
operators V which clearly have operator norm = 1. 

As for the symmetry: Clearly, —d% is symmetric. As for the other terms, 
we use the properties of the projections V ± and the operators Q± to show 
symmetry. In the case A = 0, an example of one of the terms is 



' (v 2 h)w dv k dx 



f f ^-V ± {v 2 h)V ± (v 2 k)w dv dx 



which is symmetric. (Note that in this calculation we multiply and divide 
by w since the projection operators are defined in In the case A > 0, 
we use the 'almost-symmetry' property of Q± proved in Lemma [23H5]) : 
(Q± m i n ) w = ( m > Q±™) w i which is actually a symmetry property if m 
and n are functions of x alone. This is precisely the case here. 

The discreteness of the spectrum of A^, i = 1,2, A > 0, is due to the fact 
that —d 2 defined in L p has discrete spectrum, and all other terms are 
bounded symmetric operators in L p , and, therefore, are relatively compact 
perturbations of —d\. Thus, according to the Kato-Rellich theorem 
V, §4.1], .4^ and A 2 , A > 0, are selfadjoint with domains H PQ and Hp, 
respectively. Since —d%. has pure point spectrum, so do the operators Af, 
in view of Weyl's theorem [26j IV, Theorem 5.35]. 
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2. We calculate 



\Aih-Alh\\ Ll = 



<C||Qi/ l -P ± ^|| w ^O, 



as A — > 0, by Lemma I2~5t f2"j) . and using the fact that |jujf | are bounded by 
w. Here, C = J2± sn Px (I \^t\ dv). The proof for A% follows in precisely 
the same way. 



3. Let h(x) G Hp. Then one has 



\A*h-Alh\\ Ll = 



J ^{Q^h-Qlh) dv 



< Cj2\\Q± h -Q± h \U < C|lnA-lna| \\h\ 



Again, the proof for A% follows in precisely the same way. 



4. 



\A\h+dlh\\ Ll 



<C^2\\Qih-h\\ w ^0 



as A — > oo, by Lemma |231| 4")) . Here C = J2± SU P X (J I Me I 

5. This is clearly true by verifying that A®1 = A$l = since Q A 1 = 1. 

6. The proof for A% is straightforward: By Lemma [2.5t [Tj). Q± both have 
operator norm 1. Thus, the only unbounded terms in A$ arc —d 2 and A 2 , 
both positive, and the claim follows. 

The proof for A\ is more delicate. Since the following calculations are 
quite lengthy, we drop the ±, but the same proof still stands when both 
species are considered. We need to show that for A sufficiently large and 
I 2 



for he Hi 



< (A$h, h) = [ {d x hf dx + [ [ /j, e (Q x h - h) h dv dx. (7.1) 
Jo Jo Jr 2 

Since the Laplacian on a periodic domain has a certain spectral gap G = 
G{P) > 0, we need to show that 

/ / Me (Q X h -h)hdv dx < G. (7.2) 
Jo Jr 2 

Letting — oo < /3 < be some constant to be chosen later, we may write 
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the left-hand-side of (|7.2[) as 



nMe (Q X h -h)h dv dx 
j 2 



< 



Ae As / He {h(X(s, x, v) — h{x)) h(x) dv dx ds 

-oo JO JR 2 



(7.3) 



,P , 

Ae As / He (h(X(s, x, v) — h{x)) h(x) dv dx ds 
p Jo Jm 2 



which we denote by / and II respectively. We first analyze the term I: 
\I\ = 



Ae As / / fj, e (h(X (s , x , v) — h(x)) h(x) dv dx ds 



< 







Ae As ds sup 



e~ x ^ sup 

-oo<s<^ 



'0 JS. 2 



[i e (h(X(s,x,v)) — h(x))h(x) dv dx 
c^/i(£) d£ h(x) dv dx 



Now, given s < define 

s := max {t \ X(t, x, v) £ {X(s, x, v) + mP, m £ Z}} . (7.4) 

Then due to the periodicity in the x variable, we can replace s by s in the 
last integral, and we therefore have: 



|/| < e~ X ^ sup 

— ao<s<P 



P r ( rX{s,x,v) 

He 



daHQ d u M*) dvdx 



<e~ X ^\\h\\ L 2 SUp 

-oo<s</3 \ Jo 



dv 



(i + (v)Y 



2 \ V2 

dx 



Since \X\ = \Vi\ < 1 we further have that \X(s,x,v) — x\ < P. We can 
therefore finally estimate 



m<e- A ^||/.|| i2 sup 

-oc<s<l3 \ Jo 



pX(s,x,v) \ 



1/2 



(i + (v)Y 



dx 



I L 2 SUp 

-oo<s</3 



V 



a 


/ pi/2 (7 




'K 2 \Jo 



1/2 



dv 



(i + (v)y 



\ V2 

dx 



p e - A ^||^!| L2 ||^|| L2 ||(i + («))- Q | 



L 1 • 
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Now we turn to estimating the term II: 



1^1 



Ae 



As 



< 



o Jr 2 



Ii e {h{X(s, x, v) — h(x)) h(x) dv dx ds 



Ae ds ) sup 

/3<s<0 



dv 



\h{X(s,x,v) - h{x)\ \h{x)\ d 



< WHl* sup 

/3<s<0 Wo 



< II^-Hl 2 sup 

P<s<0 \ Jo 



\h(X(s,x,v) - h(x)\ 

\o e h(o\ da 



dv 



{l + (v))° 
dv 

(i + (v)r 



2 \ 1/2 

dx 

2 \ V2 

dx 



We now again use the fact that \X\ = \Vi\ < 1, deducing that \X(s, x, v) 
x\ < \ s \ < Therefore if L8| < P, we have 



II < \\h\\i,2 sup 

f3<s<Q \ Jo 



< sup 

P<s<0 



K 2 



\ 1/2 



,1/2 



X(s,a;,i>) 



1/2 



\dM)\ ^ 



dv 



(i + (v)y 



2 

dx 



< \\h\\ L -2 sup 

P<s<0 



\ 



rP \ 1/2 

\s\ 1/2 [ I iw)i 2 da) 



dv 



(1 + (v))° 




= p 1 / 2 l|/ l || i2 ||a^|| L2 ||(i + ( v ))- Q || il |/3| 1 / 2 . 

Combining our estimates for the two terms / and II we have: 
p 



f f Ve{Q X h-h)hdvdx <\\(l + (v))- a \\ L1 \\h]\ L 4d x h\\^(Pe- m +P 1/2 W\ 1/2 ) 

< ||(1 + {v))- a \\ L iKG (Pe~ x W + P 1 / 2 ^ 1 / 2 ) 

where K is the best constant given by Poincare's inequality on [0,P]. 
Recalling (|7.2p . we need to choose f3 and A such that 



pl/2 e - m + i^jl/2 < 



(7.5) 



||(l + («))-«|| Ll /vPl/2- 

This is easily satisfiable by letting |/3| = A -1 ' 2 and taking A sufficiently 
large that only depends on P and a. We conclude that for for any < 
7 < G, there exists A > such that for all A > A, A* > 7 > 0. 
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□ 

Lemma 7.2 (Properties of B X ,C X ,V X ). Let < A < oo. 

1. B x maps L 2 P — > L 2 P with operator bound independent of A. 

2. For all h(x) E L%, ||S A /i|| L 2 p -> as A -> 0. The same is true for C A ,X> A . 

3. If a > 0, then ||S A -S t7 || = 0(\X — a\) as A -> cr, where || • || is the operator 
norm from Lp to Lp. The same is true for C A ,P A . 

4. For all ft(a;) £ L 2 P , \\B x h\\ L 2 -> as A -> oo. The same is true for C X ,V X . 
Proof. 1. Let h,k E L 2 P . 



fi p dv I h{x)k(x) dx 



E 



fJ,fQ±(v2h) cfo &(a;) dx 



/ + // 



where 



/Up dv I h(x)k{x) 



dx<^sup(| | M ±| dv) ||ft|| L , ||ft|| L , <C||%,||fc|| £ , 



and 



77 - Ej^d^Si^^)^) |^ 



< 



< 



1/2 



|fc(x)r dxj 

1/2 



y o J <*« (| 1^1 |si(^)f d«) 

Wsup/| 



2. To show that B x -> strongly in Lp we use the fact that -> 7 ?± 
strongly as shown in Lemma 12.51 We consider the two terms that make 
up B x separately. As for the first term of B x , J2± f ^p{e,p ± ) dv: Since 
(e,P + ) = -fi~(e,-p + ) one has 



[Mp (e,P + ) + M p (e,p )}dv= / ((u) , -v 2 - + fj, p {(v) ,v 2 - ip )] dv. 



The change of variables v 2 — > —v 2 applied to y, v (e, — p + ) yields cancella- 
tion. Thus the first term of B x is 0. (We note that even though this term, 
which does not depend upon A, vanishes, we keep it in this paper, as it 
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arrises naturally from Maxwell's Equations.) The second term is slightly 
more involved. As mentioned above, since Q± — > V strongly in 1/^, we 
have that 

± ± ■> 

strongly in Lp as A — > 0. Therefore, we now show that J2± J ^f'P ± {v2h)dv = 
0. Observe that since E° = 0, g(x, v) £ ker D~ if and only if g(x, —v) £ 
kerD + . Therefore V [g(x,v)] = g(x,v) if and only if V + [g(x, — v)] = 
g(x,—v). We conclude that V~ [g(x, v)](x, —v) — g(x, — v) = V + [g(x, — v)](x, v). 
Hence the second term of B x h converges to 



[nfV + [v2h](x, v) + fi e V [v2h](x,v)) dv = j {ji^V [-v 2 h](x, —v) + ^, e V [v 2 h](x, v)) dv 

(— Me [7)2/1] (x, —v) + \J^,V~ [v2h](x, v)) dv 



where in the second equality we used the linearity of V~ . Now, wc know 
that /Lt+(e,p + ) = d e fi + ((v) , U2 + V' ) = ®eH~({v) ,~v 2 — ip°). Therefore, by 
changing variables v —> — v in the first integrand we get exact cancellation 
once again. Therefore, indeed, 23 A — s- strongly in Lp. 

The fact that ||C A (6)|| L 2 -> (and similarly for V x ) holds since (i) as A ^ 
0, || Q±m— V ± m\\ w — > 0, as we have seen in LemmaES] (ii) fi is even in vi, 
and, finally, (iii) the projection operators V ± preserve parity with respect 
to the variable Vi by Lemma 12^1 Thus ^ ± J fJ^V (vi) dv = 0. Combin- 
ing these three facts, we get that ||C A (b) \\ L 2 \\bJ2± J Me ^ ?± ( 1 ^i) ^lli 2 = 
0. 

3. The proof of this fact is identical to the proof for A x . 

4. By the triangle inequality: 



\B X h\\ T ,2 < 



h J2f + <W0 dV + E / £ (Q±& h ) - *2h) dv 

J- J J- J 



I+II 



where / vanishes due to (|2.13[) . The term 77 is controlled as follows: 



11 



Me {Qi(v2h)-v 2 h)\dv\ dx 



1/2 



< 



< 



{/ f \^\\Q±^2h)-v 2 h\ 2 dv^dx 



1/2 



E 



sup y \fif\dv |y |//^ \Qi(v 2 h) - v 2 h\ 2 dv^dx^ <Cj2\\Q±(^ h )~^ h \\ w ^° 
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as A — > oo by Lemma 12.51 



As for C x , we summarize the argument: 



53 f dv 


-> 


X] / A^e "1 ^ 


= 


± J 




± ^ 


i 2 p 



as A — > oo by Lemma 12.51 and the fact that [i is even in v± . The same 
proof holds for T> x . 



Lemma 7.3 (Properties of l x ). Let < A < oo. 

1. l x -^/°asA^0. 

2. l x is uniformly bounded in A. 

Proof. 1. This is an immediate consequence of the properties of Q±. As we 
have demonstrated how to use these properties above, we do not repeat 
the proof. 

2. This is due to the fact that Q± have operator norm = 1, the integrability 
properties of fif, and the fact that |£>i| < 1. 



The following lemma lists the important properties of M x - all of which are 
inherited directly from the properties of the various operators it is made up of, 
as listed in Lemma I7TT1 

Lemma 7.4 (Properties of M x ). To simplify notation, we write u T for a generic 
element (0, if), b) £ Hp x Hp x E. 

1. For all A > 0, M x is selfadjoint on Lip x Lp x K with domain Hp x Hp x R. 

2. For all u T £ Hp x Hp x R, ||A4 A u - 7W°u|| L 2 pXL 2 jXi 2 p -> as A -> 0. 

3. If <7 > 0, then ||A1 A — CT 1 1 — > as A — > tr, where || • || is the operator 
norm from Hp x Hp x R to Lp x Lp x Lp. 

Proof. 1. This is true due to the structure of Ai x and the selfadjointness of 
A x and A\ for A > 0. 

2. We note that 



□ 



□ 
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Thus 



\\(M x -M Q )u\\ L2pXL%xL% < ||(-^+^)0!| + ||6Vll + l|C A 6|| 

+ \\{B x )*cf ) \\ + \\{A x -A 2 )i>\\ + \\V x b\\ 

+ || (C X Y <j>\\ + || {V x )*^\\ + || (-P (A 2 - l x ) - Pl°) fell . 

where all norms on the right hand side are in L 2 P . Since each of the terms 
on the right hand side tends to as A — > 0, we are done. 

3. We want to show that j| A4 X — A4 a \\ — >• as A — s- a. This, again, is true by 
virtue of the fact that this is true for each of the entries of M. x separately. 

□ 

Acknowledgements: This paper is part of the author's Ph.D. thesis, super- 
vised by Professor Walter Strauss. It is a great pleasure to thank Professor 
Strauss for his patient and kind guidance, his generous advice and the numer- 
ous hours he spent going over all details during the writing of this work. 



References 

[1] Bernstein, I. B.; Frieman, E. A.; Kruskal, M. D.; Kulsrud, R. M., An energy 
principle for hydromagnetic stability problems. Proc. Roy. Soc. London. Ser. 
A. 244 (1958), 17-40. 

[2] Bernstein, I., Greene, J., Kruskal, M., Exact nonlinear plasma oscillations. 
Phys. Rev. 108, 3, 17-40, (1958). 

[3] Courant, Richard and Hilbert, David, Methods of Mathematical Physics, 
Vol. 1, Interscience Publishers 1966. 

[4] Davidson, Ronald C, Physics of Nonneutral Plasmas, Addison- Wesley 
1990. 

[5] Davidson, Ronald C. and Hong, Qin, Physics of intense charged particle 
beams in high energy accelerators, World Scientific, 2001. 

[6] Dendy, R. O., Plasma Dynamics, Oxford Science Publications, 1990. 

[7] DiPcrna, Ronald, and Lions, Pierre-Louis, On the Cauchy Problem for 
Boltzmann Equations: Global Existence and Weak Stability, Ann. Math. 
130 (1989), 321-366. 

[8] DiPerna, Ronald, and Lions, Pierre-Louis, Global Weak Solutions of 
Vlasov-Maxwell Systems, Comm. Pure Appl. Math., 42 (1989), 729-757. 

[9] Friedberg, J. P., Ideal Magnetohydrodynamics, Plenum Press, 1987. 



43 



[10] Gardner, C.S. Bound on the energy available from a plasma, Phys. Fluids, 
6 (1963), 839-840. 

[11] Glassey, Robert, The Cauchy Problem in Kinetic Theory, SIAM, 1996. 

[12] Glassey, Robert and Strauss, Walter, Singularity formation in a collision- 
less plasma could occur only at high velocities, Arch. Ration. Mech. Anal. 
92 (1986), no. 1, 59-90. 

[13] Glassey, Robert and Schaeffer, Jack, Global existence for the relativistic 
Vlasov- Maxwell system with nearly neutral initial data, Comm. Math. Phys. 
199 (1988), no. 3, 353-384. 

[14] Glassey, Robert and Schaeffer, Jack, On the "one and one-half dimen- 
sional" relativistic Vlasov- Maxwell system, Math. Methods Appl. Sci. 13 
(1990), no. 2, 169-179. 

[15] Glassey, Robert and Schaeffer, Jack, The "two and one-half dimensional" 
relativistic Vlasov-Maxwell system, Comm. Math. Phys. 185 (1997), 257- 
284. 

[16] Goedbloed, Hans and Poedts, Stefan, Principles of Magnetohydrodynamics 
: With Applications to Laboratory and Astrophysical Plasmas, Cambridge 
University Press, 2004. 

[17] Grad, Harold, The guiding center plasma. 1967 Proc. Sympos. Appl. Math., 
Vol. XVIII pp. 162-248, Amcr. Math. Soc, Providence, R.I. 

[18] Guo, Yan, Stable magnetic eguilibria in collisionless plasma, Comm. Pure 
Appl. Math., Vol L, 0891-0933 (1997). 

[19] Guo, Yan, Stable magnetic equilibria in a symmetric plasma, Commun. 
Math. Phys., 200, 211-247 (1999). 

[20] Guo, Yan and Strauss, Walter, Instability of periodic BGK equilibria, 
Comm. Pure Appl. Math. Vol XLVIII, 861-894 (1995). 

[21] Guo, Yan and Strauss, Walter, Unstable oscillatory-tail waves in collision- 
less plasmas, SIAM J. Math. Anal., 30, no. 5, 1076-1114 (1999). 

[22] Guo, Yan and Strauss, Walter, Unstable relativistic BGK waves, Comput. 
and Appl. Math., 18, no. 1, 87-122 (1999) 

[23] Guo, Yan and Strauss, Walter, Unstable BGK solitary waves and collision- 
less shocks, Commun. Math. Phys. 195, 249-265 (1998). 

[24] Guo, Yan and Strauss, Walter, Magnetically created instability in a colli- 
sionless plasma, J. Math. Pures. Appl., 79, no. 10, 975-1009 (2000). 

[25] Holm, D., Marsden, J., Ratiu, T., Weinstein, A., Nonlinear stability of fluid 
and plasma equilibria, Physics Reports, 123, no. 1-2, 1-116 (1985). 



44 



[26] Kato, Tosio, Perturbation Theory for linear operators , (2nd edition) 
Springer 1976. 

[27] Kruskal, Martin, Hydromagnetics and the theory of plasma in a strong 
magnetic field, and the energy principles for equilibrium and for stability. 
1960 La theorie des gaz neutres et ionises (Grenoble, 1959) pp. 251-274 
Hermann, Paris; Wiley, New York 

[28] Kulsrud, Russell, General stability theory in plasma physics. 1964 Advanced 
Plasma Theory (Proc. Intcrnat. School of Physics "Enrico Fermi", Course 
XXV, Varenna, pp. 54-96 Academic Press, New York 

[29] Langmuir, Irving, Oscillations in Ionized Gases, Proc. Natl. Acad. Sci.. Vol 
14, 8 627-637 (1928). 

[30] Lavel, G., Mercicr, C. and Pcllat, R.M., Necessity of the energy principle 
for magnetostatic stability, Nuclear Fusion 5 , 156-158 (1965). 

[31] Lin, Zhiwu, Instability of periodic BGK waves, Math. Res. Letts., 8, 521- 
534 (2001). 

[32] Lin, Zhiwu, Nonlinear instability of periodic waves for Vlasov-Poisson sys- 
tem, Comm. Pure Appl. Math. Vol LVIII, 505-528 (2005). 

[33] Lin, Zhiwu, and Strauss, Walter, A sharp stability criterion for the Vlasov- 
Maxwell system, Invent. Math. 173, 497-546 (2008) 

[34] Lin, Zhiwu, and Strauss, Walter, Linear stability and instability of relativis- 
ts Vlasov-Maxwell systems Comm. Pure Appl. Math. 60, 724-787 (2007) 

[35] Lin, Zhiwu and Strauss, Walter, Nonlinear stability and instability of rela- 
tivistic Vlasov-Maxwell systems, Comm. Pure Appl. Math. 60, 789-837 
(2007) 

[36] Newcomb, W. A., Lagrangian stability of MHD fluids, Nucl. Fusion. Suppl. 
2, 451. 

[37] Nicholson. D. R., Introduction to Plasma Theory, Wiley, 1983. 

[38] Parks, George, Physics of space plasmas, Second Edition, Westview Press, 
2004. 

[39] Penrose, Oliver, Electrostatic Instabilities of a Uniform Non-Maxwellian 
Plasma , Phys. Fluids, 3, 258, 258-265 (1960). 

[40] Pfaffclmoscr, K., Global classical solutions of the Vlasov-Poisson system in 
three dimensions for general initial data, J. Diff. Eqns. 95, 281-303 (1992). 

[41] Sturrock, P., Plasma Physics - An introduction to the theory of astrophysi- 
cal, geophysical, and laboratory plasmas. Cambridge University Press, 1994. 



45 



[42] Taylor, J.B., Plasma containment and stability theory, Proceedings of the 
Royal Society of London. Series A, Mathematical and Physical Sciences, 
Vol. 304, No. 1478. (1968), pp. 335-360. 

[43] Trivelpiece, A. W. and Krall, N. A., Principles of Plasma Physics, 
McGraw-Hill, 1973. 

[44] Vlasov, Anatoly, On Vibration Properties of Electron Gas (Russian), J. 
Exp. Theor. Phys., Vol. 8, No. 3. (1938), p. 291. 



46 



